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Abstract —A semideiinite programming (SDP) relaxation glob¬ 
ally solves many optimal power flow (OPF) problems. For other 
OPF problems where the SDP relaxation only provides a lower 
bound on the objective value rather than the globally optimal 
decision variables, recent literature has proposed a penalization 
approach to find feasible points that are often nearly globally 
optimal. A disadvantage of this penalization approach is the 
need to specify penalty parameters. This paper presents an alter¬ 
native approach that algorithmically determines a penalization 
appropriate for many OPF problems. The proposed approach 
constrains the generation cost to be close to the lower bound 
from the SDP relaxation. The objective function is specified 
using iteratively determined weights for a Laplacian matrix. This 
approach yields feasible points to the OPF problem that are 
guaranteed to have objective values near the global optimum 
due to the constraint on generation cost. The proposed approach 
is demonstrated on both small OPF problems and a variety of 
large test cases representing portions of European power systems. 

Index Terms —Optimal power flow, Semideiinite optimization. 
Global solution 

I. Introduction 

HE optimal power flow (OPF) problem determines an 
optimal operating point for an electric power system 
in terms of a specified objective function (typically gener¬ 
ation cost per unit time). Equality constraints for the OPF 
problem are dictated by the network physics (i.e., the power 
flow equations) and inequality constraints are determined by 
engineering limits (e.g., voltage magnitudes, line flows, and 
generator outputs). 

The OPF problem is non-convex due to the non-linear power 
flow equations, may have local optima Cl, and is generally 
NP-Hard El, even for relatively simple cases such as tree- 
topologies El. There is a large literature on solving OPF 
problems using local optimization techniques (e.g., successive 
quadratic programs, Lagrangian relaxation, heuristic optimiza¬ 
tion, and interior point methods a, 0). These techniques 
are generally well suited to solving large problems. However, 
while local solution techniques often find global solutions ||6|, 
they may fail to converge or converge to a local optimum m, 
Q. Furthermore, they are unable to quantify solution optimal¬ 
ity relative to the global solution. 

There has been significant recent research focused on con¬ 
vex relaxations of OPF problems. Convex relaxations lower 
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bound the optimal objective value and can certify infeasibility 
of OPF problems. These relaxations also yield the globally 
optimal decision variables for many OPF problems (i.e., the 
relaxations are often “exact”). Second-order cone program¬ 
ming (SOCP) relaxations can globally solve OPF problems 
for radial networks that satisfy certain non-trivial technical 
conditions 0. Semidefinite programming (SDP) relaxations 
globally solve a broader class of OPF problems El, lH, iflOl . 

Recently, the SDP relaxation has been generalized to a fam¬ 
ily of “moment” relaxations using the Lasserre hierarchy for 
polynomial optimization lfTT]| - lfT3ll . With increasing relaxation 
order, the moment relaxations globally solve a broader class 
of OPF problems at the computational cost of larger SDPs. 

By exploiting network sparsity and selectively applying 
the computationally intensive higher-order constraints, the 
moment relaxations are capable of globally solving larger OPF 
problems m, including problems with several thousand buses 
representing portions of European power systems Ca. On¬ 
going efforts include further increasing computational speed. 
Recent work includes implementing the higher-order moment 
constraints with a faster SOCP formulation ifTbll and develop¬ 
ment of a complex version of the Lasserre hierarchy ini. 

As an alternative to the moment relaxations, other literature 
has proposed an objective function penalization approach for 
finding feasible points that are near the global optimum for the 
OPF problem lITSl . 1191 . The penalization approach has the 
advantage of not using potentially computationally expensive 
higher-order moment constraints, but has the disadvantage of 
requiring the choice of appropriate penalization parameters. 
This choice involves a compromise, as the parameters must 
induce a feasible solution to the original problem while 
avoiding large modifications to the problem that would cause 
unacceptable deviation from the global optimum. 

The penalization formulation in the existing literature lIT^ 
generally requires specifying penalty parameters for both the 
total reactive power injection and apparent power flows on 
certain lines. Penalty parameters in the literature range over 
several orders of magnitude for various test cases, and existing 
literature largely lacks systematic algorithms for determining 
appropriate parameter values. Recent work ifTsll proposes a 
“momentH-penalization” approach that eliminates the need to 
choose apparent power flow penalization parameters, but still 
requires selection of a penalty parameter associated with the 
total reactive power injection. 

This paper presents an iterative algorithm that builds an 
objective function intended to yield near-globally-optimal so¬ 
lutions to OPF problems. The algorithm is applicable for 
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cases where the SDP relaxation is not exact but has a small 
“relaxation gap” (i.e., the solution to the SDP relaxation has 
an objective value that is close to the true globally optimal 
objective value). The proposed algorithm first solves the SDP 
relaxation to obtain a lower bound on the optimal objective 
value. For many practical OPF problems, this lower bound is 
often very close to the global optimum. The proposed approach 
modifies the SDP relaxation by adding a constraint that the 
generation cost must be within a small percentage (e.g., 0.5%) 
of this lower bound. This percentage is the single externally 
specified parameter in the proposed approach. 

This constraint on the generation cost provides freedom to 
specify an objective function that aims to obtain a feasible 
rather than minimum-cost solution for the OPF problem. In 
other words, we desire an objective function such that the 
SDP relaxation yields a feasible solution to the original non- 
convex OPF problem, with near-global optimality ensured by 
the constraint on generation cost. 

This paper proposes an algorithm for calculating an appro¬ 
priate objective function defined using a weighted Laplacian 
matrix. The weights are determined iteratively based on the 
mismatch between the solution to the relaxation and the power 
flows resulting from a related set of voltages. The paper will 
formalize these concepts and demonstrate that this approach 
results in near global solutions to many OPF problems, in¬ 
cluding large test cases. Like many penalization/regularization 
techniques flSl . lIT^ . the prmosed approach is not guaranteed 
to yield a feasible solutionis As supported by the results for 
several large-scale, realistic test cases, the proposed algorithm 
broadens the applicability of the SDP relaxation to achieve 
operating points for many OPF problems that are within 
specified tolerances for both constraint feasibility and global 
optimality. 

There is related work that chooses the objective function of 
a relaxation for the purpose of obtaining a feasible solution for 
the original non-convex problem. For instance, ll20l specifies 
objective functions that are linear combinations of squared 
voltage magnitudes in order to find multiple solutions to the 
power flow equations. Additionally, ED proposes a method 
for determining an objective function that yields solutions to 
the power flow equations for a variety of parameter choices. 
The objective function in ||2TI is defined by a matrix with three 
properties: positive semidefiniteness, a simple eigenvalue of 0, 
and null space containing the all-ones vector. We note that the 
weighted Laplacian objective function developed in this paper 
is a special case of an objective function that also has these 
three properties. 

This paper is organized as follows. Section [II] introduces 
the OPF formulation studied in this paper. Section Hill reviews 
the SDP relaxation from previous literature. Section |IV] de¬ 
scribes the Laplacian objective function approach that is the 
main contribution of this paper. Section |V] demonstrates the 
effectiveness of the proposed approach through application 

^ The moment+penalization method in [B] applies the Lasserre hierarchy 
for polynomial optimization to OPF problems that are augmented with a 
rea ctiv e power penalty. For sufficiently high relaxation orders, the approach 
in QU is guaranteed to yield a feasible solution. However, solving relaxations 
from the Lasserre hierarchy can be computationally challenging. 


to a variety of small OPF problems as well as several large 
test cases representing portions of European power systems. 
Section IVTl concludes the paper. 

II. Optimal Power Flow Problem 

We first present an OPF formulation in terms of complex 
voltage coordinates, active and reactive power injections, and 
apparent power line flow limits. Consider an n-bus system 
with ni lines, where Af = { 1 ,..., n} is the set of buses, Q 
is the set of generator buses, and £ is the set of lines. The 
network admittance matrix is Y = G + jB, where j denotes 
the imaginary unit. Let P^k +jQDk represent the active and 
reactive load demand and 14 = Vdk +jVqk the voltage phasors 
at each bus k G Af. Superscripts “max” and “min” denote 
specified upper and lower limits. Buses without generators 
have maximum and minimum generation set to zero. Let Ck 2 , 
Cfci, and Cko denote the coefficients of a convex quadratic cost 
function for each generator k G Q. 

The power flow equations describe the network physics: 

n 

Pck —^dk ^ ^ ^ikf^qi) 

n 

+ Vqk (BifcVdi -|- GikVqi) -\- Pok (1^) 

i=l 
n 

Qck —^dk ^ ^ ( ^ik^di Gik^qi) 
n 

+ Vqk {GikVdi — 'RikVqi) Qok- (lb) 

i=l 

We use a line model with an ideal transformer that has a 
specified turns ratio : 1 in series with a II circuit with 

series impedance Rim +jXim (equivalent to an admittance of 
gim +jbim = shunt susceptance ibsh,im- 

The line flow equations are 

Pl-m = {Vil -G Vql) gimlnm 

(VdtVdm P VqtVqm) iflm Sin {Otm} Qlm COS l^lm 

-p 14^ P^dm) (, 91 m, SlU (Olm) P blm COS (Him}} j '^Im 

(2a) 

Pml — ifPdm P ^qm) Qlm 

(14/ P^dm P Pqt Pqm} (Qlm COS (Olm} P blm SIH (Olm}} / Rm 
P (PdlPqm PqlPdm} (Qtm slu (Olm} (//m COS (Olm}} IRm 

(2b) 

Qlm = - {Vil P Vq]} (^blm P /tL 

“h (VdiVdm V^iVqm) (^*^771 COS (^Zm) H“ Qlm Sill j'^lm 

“h (y^dl'^qm ^ql^dm) {Qlm COS (^^771) Sill (^irri)) l'l~lm 

(2c) 

Qml = - {VL P {him P 

~\~ {Vdl^dm ~\~ V^TTi) (^^TTi COS {Olm) Qlm sill j'^lm 

-|- VqlVdm^ {Qlm COS ^Im sill j'l'lm- 

(2d) 

The classical OPF problem is then 
min Ck 2 PGk + CkiPok P Cfco subject to ( 3 a) 

P^t<PGk<P^r ^kGff ( 3 b) 
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Qsr < Qck < Qgr 

ykeAf 

(3c) 

iVk^n^<vi + Vqi<iVk^n^ 

WkGAf 

(3d) 

{Plmf + {Qlmf < {SZn" 

y{l,m) e £ 

(3e) 

{Pmlf + {Qmlf < 

V {1, m) £ C 

(3f) 

II 

O 


(3g) 


Constraint (3g i sets the reference bus angle to zero. 


III. Semidefinite Relaxation of the OFF Problem 


This section describes an SDP relaxation of the OPF 
problem adopted from ||3, IfTOl . Il22l . Let denote the 
standard basis vector in K". Dehne Yk = ekeJ.Y, where (•)^ 
indicates the transpose operator. 

Matrices employed in the bus power injection, voltage 
magnitude, and angle reference constraints are 


Yk 

Yk 

Mk 

Nfc 


Re{Yk + Y^) 
lm{Yk-Y^) 
lm{Yk+Y^) 
Re{Yj-Yk) 

0 

0 Ckel 




ImiY^-Yk) 
Re{Yk +Y^)_ 
Re{Yk-Y^) 
ImiYk + Y^) 


0 0 

0 Ckcl 


(4a) 

(4b) 

(4c) 

(4d) 


where Re( ) and Im( ) return the real and imaginary parts, 
respectively, of a complex argument. 

Define fi as the standard basis vector in and define: 


Qm — iSlm cos {Olm) sin (5a) 

^mi ~ {,9im COS {OifYi') -\- bi^ sin / (2 t/^) (5b) 

^Im — idlm sin {Olm) “t“ COS (5c) 

^ml — {dim sin {Olm') COS / {^'^Im) • (5d) 

Matrices employed in the line flow constraints can then be 
written: 


Z/m. — 


9ln 


' Im 


[flfl + //+„/;+„) 


- Clra {fJL + fmfl + fl+nfli+n + fm+nfj+n) 


Sim 


{flfl+n + fm+nfl “ fl+nfl “ fm9+n) 


( 6 a) 


dim i^fmfm fm+nfm+n) 

- Cml {flfl + fmfl + fl+nfm+n + fm+nfl+n) 

- Sml (fl+nfl + /m/z+n “ flfl+n “ fm+nfl) (6b) 

^bim bsh,lr, 


— 


2 t 2 


{flfl + fl+nfl+n) 


+ Clm {fifl+n + fm+nfl “ fl+nfl “ fmfl+n) 

+ Sim {flfl + fmfl + fl+nfm+n + fm+nfl+n) (6c) 

Zml = -( blm + ^ 


{fmfl + fm +nf m+n) 

+ C^Z {fl+nfl + fmfl+n “ flfl+n “ fm+nfl) 
Sml {flfm~^fmfI fl+nfm+n fm+nfl+n) ' 


( 6 d) 


Define the vector of voltage components: 

X=[Vdl Vd2 ... Vdn Vgi Vg2 ... VgnY (7) 
and the rank-one matrix: 

W = a;xL ( 8 ) 

The active and reactive power injections at bus k are 
tr (YfcW) and tr (Y^W), respectively, where tr (.) indicates 
the matrix trace operator. The square of the voltage magnitude 
at bus k is tr (M^W). The constraint tr (NiW) = 0 sets the 
reference angle. 

Replacing the rank-one requirement from ([ 8 ]l by the less 
stringent constraint W ^ 0, where ^ 0 indicates positive 
semidehniteness, yields the SDP relaxation of ([^i: 


mm > ak subject to 
'W,a,Pc 

keg 

Pck = tr (Yk'W) + Pok 

Pat < Pgu < P^r 


Qak < tr (Y,W) + Quk < QgT 

(pr") 

<tr(M;,W) < {yr") 

tr(NiW) = 0 


(1 — CklPGk — CkO + O-k) 

> 

,1 + CklPGk + CkO — Cik) 

2yJCk2PGk \ 

omax \ 

l^lm 
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^Itn ^ 

[tr (Zi^W). 
rtr(ZwW)- 

2 

[tr 

2 

W ^ 0. 




Vfe G Af 
Vfc € Af 
Vfc € Af 

Vfc € Af 


Mkeg 

V {I, m) € C 

V {I, m) £ C 


(9a) 

(9b) 

(9c) 

(9d) 

(9e) 

(9f) 

(9g) 

(9h) 

(9i) 

(9j) 


The generation cost constraint is implemented using 
the auxiliary variable ak and the SOCP formulation in (9gl. 
The apparent power line flow constraints (1^ and (O are 
implemented with the SOCP formulations in ( l9h] i and ( l9i] l. 
See ifTOll for a more general formulation of the SDP relaxation 
that considers the possibilities of multiple generators per bus 
and convex piecewise-linear generation costs. 

Note that rather than explicitly constraining the reference 
angle, the constraint ( [3^ can be used to eliminate the variable 
Vqi from the problem. Eliminating Vqi removes the n-l-1 row 
and column from W, with corresponding modihcations to all 
matrices in and removal of ( l9fl ). This approach is often 
numerically superior to explicitly constraining the reference 
angle as in ( l9ft . 

If the condition rank (W) = 1 is satished, the relaxation 
is “exact” and the global solution to Q is recovered using 
an eigen-decomposition. Let A be the non-zero eigenvalue 
of a rank-one solution W to (|9ll with associated unit-length 
eigenvector 77 . The globally optimal voltage phasor is 


E* = \/A {r]i.,n + j??(„+l):2n) (10) 


where subscripts denote vector entries in MATLAB notation. 

The computational bottleneck of the SDP relaxation is 
the constraint < 11 ^ which enforces positive semidehniteness 
for a 2n X 2n matrix. Solving the SDP relaxation of large 
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OPF problems requires exploiting network sparsity. A matrix 
completion decomposition exploits sparsity by converting the 
positive semidefinite constraint on the large W matrix (|^ to 
positive semidefinite constraints on many smaller submatrices 
of W. These submatrices are defined using the cliques (i.e., 
completely connected subgraphs) of a chordal extension of the 
power system network graph. See IfTOll . Il22l . ||2^ for a full 
description of a formulation that enables solution of (|9ll for 
systems with thousands of buses. 

IV. Laplacian Objective Function 

The SDP relaxation in Section|ni]globally solves many OPF 
problems ||2|, ifTOll . However, there are example problems for 
which the SDP relaxation fails to yield the globally optimal 
decision variables (i.e., the solution to the SDP relaxation does 
not satisfy the rank condition ([8]l). This section proposes an 
approach for finding feasible points near the global optima 
of many problems for which the lower bounds from the SDP 
relaxation are close to the globally optimal objective values. 

The proposed approach constrains the generation cost to be 
close to the lower bound obtained from the SDP relaxation. 
This enables the specification of an objective function based 
on a weighted Laplacian matrix that yields feasible (i.e., rank- 
one) solutions to many OPF problems. An iterative algorithm 
based on line flow mismatches is used to determine the weights 
for the Laplacian matrix. 

A. Generation Cost Constraint 

The proposed approach exploits the empirical observation 
that the SDP relaxation provides a very close lower bound on 
the optimal objective value of many typical OPF problems 
(i.e., there is a very small relaxation gap). For instance, 
the SDP relaxation gaps for the large-scale Polish UtM and 
PEGASE ll^ systems, which represent portions of European 
power systems, are all less than 0.3%lxl Eurther, the SDP 
relaxation is exact (i.e., zero relaxation gap) for the IEEE 
14-, 30-, 39-, 57-bus systems, the 118-bus system modified 
to enforce a small minimum line resistance El, and several 
of the large-scale Polish test cases 1110 (See m, m for 
case descriptions.) Numerical experiments also demonstrate 
that the SDP relaxation is exact for a variety of test cases with 
multiple local optima (e.g., WB2, WB3, WB5mod, and the 
22- and 30-bus loop systems in lIZTl l. To further demonstrate 
the capabilities of the SDP relaxation, 1000 modified versions 
were created for each of the IEEE 14-, 30-, 39-, and 57-bus 
systems using normal random perturbations (zero-mean, 10% 
standard deviation) of the load demands and power generation 
limits. The SDP relaxation was exact (or proved infeasibility) 

^To obtain satisfactory convergence of the SDP solver, these systems are 
pre-processed to remove low-impedance lines (i.e., lines whose impedance 
values have magnitudes less than 1 X 10“® per unit) as in OH 

^These relaxation gaps are calculated using the objective values from 
the SDP relaxation m and solutions obtained either from the second-order 
moment relaxation □3 (where possible) or from Matpower ED 

^Even the minor modifications performed when pre-processing low- 
impedance lines and enforcing minimum line resistances are not needed for 
some test cases. For instance, the SDP relaxation is exact for the Polish 
systems 2376sp, 2737sop, 2746wp, and 2746wop without modifications Qo). 


for 100% and 98.7% of the test cases derived from the 
14- and 57-bus systems, respectively. After modifications to 
enforce a 1 x 10“"^ per unit minimum line resistance, the 
SDP relaxation was exact or proved infeasibility for 81.8% 
and 81.2% of the test cases derived from the 30- and 39-bus 
systems, respectively. 

This section assumes that the lower bound provided by 
the SDP relaxation is within a given percentage 6 of the 
global optimum to the OPE problem. (Most of the examples 
in Section |V] specify <5 = 0.5%.) We constrain the generation 
cost using this assumption: 

'y ] Ck2PGk P CklPck + Cfco ^ C* (1 -|- 5) (11) 

k&g 

where c* is the lower bound on the optimal objective value 
of (H obtained from the semidefinite relaxation dUl. This 
constraint is implemented by augmenting the SDP relaxation’s 
constraints (l9b]l-([^ with 

^afc<c*(l + 5). (12) 

fcee 

If the SDP relaxation (|9l) is feasible, the feasible space 
defined by (l9b]) -([^ and (fT2li is non-empty for any choice 
of J > 00 However, if S is too small, there may not exist a 
rank-one matrix W (i.e., a feasible point for the original OPE 
problem (|2l) in the feasible space. 

The lack of a priori guarantees on the size of the relaxation 
gap is a challenge that the proposed approach shares with 
many related approaches for convex relaxations of the opti¬ 
mal power flow problem. Existing sufficient conditions that 
guarantee zero relaxation gap generally require satisfaction of 
non-trivial technical conditions and a limited set of network 
topologies 0, ESI. The SDP relaxation is, however, exact for 
a significantly broader class of OPE problems than those that 
have a priori exactness guarantees, and has a small relaxation 
gap for an even broader class of OPE problems0 

There are test cases that are specifically constructed to 
exhibit somewhat anomalous behavior in order to test the 
limits of the convex relaxations. The SDP relaxation gap 
is not small for some of these test cases. Eor instance, the 
3-bus system in ll28l . the 5-bus system in ll29l , and the 9-bus 
system in El have relaxation gaps of 20.6%, 8.9%, and 10.8%, 
respectively, and the test cases in ll^ have relaxation gaps 
as large as 52.7%. The approach proposed in this paper is 
not appropriate for such problems. Euture progress in convex 
relaxation theory is required to develop broader conditions that 
provide a priori certification that the SDP relaxation is exact 
or has a small relaxation gap. We also await the development 
of more extensive sets of OPE test cases to further explore 
the observation that many typical existing practical test cases 
have small SDP relaxation gaps. 

^Infeasibility of the SDP relxation © certifies infeasibility of the original 
OPF problem 0. 

^None of the aforementioned IEEE test cases, Polish systems, and 
PEGASE systems satisfy any known sufficient conditions for exactness of 
the SDP relaxation, but many still have zero or very small relaxation gaps. 
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B. Laplacian Objective Function 
Consider the optimization problem 
min / (W) 

w,q(,Pg (13) 

subject to da 

where / (W) is an arbitrary linear function. Any solution 
to ( foi l with rank (W) = 1 yields a feasible solution to the 
OPF problem Q within 5 of the globally optimal objective 
value due to the constraint (fa on the generation cost. This 
constraint effectively frees the choice of the function / (W) 
to obtain a feasible rather than minimum-cost solution to (l3]l. 

We therefore seek an objective function / (W) which maxi¬ 
mizes the likelihood of obtaining rank(W) = 1. This section 
describes a Laplacian form for the function / (W). Specif¬ 
ically, we consider a n; x n; diagonal matrix D containing 
weights for the network Laplacian matrix L = 
where Ainc is the ni x n incidence matrix for the network. 
The off-diagonal term Ly is equal to the negative of the sum 
of the weights for the lines connecting buses i and j, and the 
diagonal term ha is equal to the sum of the weights of the 
lines connected to bus i. The objective function is 


/ (W) = tr 


0 . 


-'nxn 

L 


w 


(14) 


The choice of an objective function based on a Laplacian 
matrix is motivated by previous literature. An existing pe¬ 
nalization approach M augments the objective function by 
adding a term that minimizes the total reactive power injection. 
This reactive power penalty can be implemented by adding the 
term 


ffetr 



W 


(15) 


to the objective function of the SDP relaxation (I9a] l. where 
£{, is a specified penalty parameter and indicates the 
complex conjugate transpose operator. In the absence of phase- 
shifting transformers (i.e., 6im = 0 \/ {l,m) G C), the matrix 

^ ^7^ is equivalent to — Im (Y) = —B, which is a weighted 
Laplacian matrix (with weights determined by the branch 
susceptance parameters bim = ) plus a diagonal 

matrix composed of shunt susceptances. 

Early work on SDP relaxations of OPF problems advo¬ 
cates enforcing a minimum resistance of for all lines in the 
network. For instance, the SDP relaxation fails to be exact for 
the IEEE 118-bus system ESj . but the relaxation is exact after 
enforcing a minimum line resistance of e^. = 1 x 10“^ per unit. 
After enforcing a minimum line resistance, the active power 
losses are given by 


tr 


Re 
— Im 



W 


(16) 


where Y^. is the network admittance matrix after enforcing 
a minimum branch resistance of e^. In the absence of phase- 


shifting transformers. 


Y,.+Yt 


is equivalent to Re (Y^), which 


is a weighted Laplacian matrix (with weights determined by 


the branch conductance parameters ) pl^^^ 

a diagonal matrix composed of shunt conductances. Since 
typical OPF problems have objective functions that increase 
with active power losses, enforcing minimum line resistances 
is similar to a weighted Laplacian penalization]^ 

The proposed objective function (fT4li is equivalent to a linear 
combination of certain components of W: 

/(W)= ^ (W„ - 

+W/+„,i+„ - 2Wz+ 

n,m-t-n T m+n^m+n) 

(17) 

where Dp m) is the diagonal element of D corresponding to 
the line from bus I to bus to. If W = xx'^ (i.e., W is a 
rank-one matrix) with x defined as in (0, then the objective 
function (fT4l l is equivalent to 

f (xx"t) — ^ ^ ^m) “t” J" 

{l,m)GC 

(18) 

Note that the Laplacian objective function (fT^ is convex 
in the voltage components 14 and Vq when the weights 
in D are non-negative. This is in contrast to the reactive 
power penalization in ifTSl : an objective function that penalizes 
reactive power injections may be non-convex in terms of the 
voltage components 14 and Vq when the network has non-zero 
shunt capacitors. 

When the SDP relaxation fails to yield the global opti¬ 
mum, the relaxation often “artificially” increases the voltage 
magnitudes to reduce active power losses. This results in 
voltage magnitudes and power injections that are feasible for 
the relaxation (|9|l but infeasible for the OPF problem (O. 
By minimizing the squared differences between the voltage 
phasors at connected buses, the Laplacian objective function 
counteracts this tendency of the SDP relaxation. Intuitively, 
the proposed approach uses the Laplacian objective function 
to balance two potentially competing tendencies: increasing 
voltage magnitudes to reduce active power losses such that 
the generation cost constraint is satisfied versus decreasing 
voltage differences to reduce the Laplacian objective function. 

From a physical perspective, the Laplacian objective’s ten¬ 
dency to reduce voltage differences is similar to both the 
reactive power penalization proposed in lITSl and the minimum 
branch resistance advocated in For typical operating con¬ 
ditions, reactive power injections are closely related to voltage 
magnitude differences, so penalizing reactive power injections 
tends to result in solutions with similar voltages. Likewise, the 
active power losses associated with line resistances increase 
with the square of the current flow through the line, which 
is determined by the voltage difference across the line. Thus, 
enforcing a minimum line resistance tends to result in solutions 
with smaller voltage differences in order to reduce losses. 

^Note that since enforcing minimum line resistances also affects the 
power injections and the line flows, the minimum line resistance cannot be 
solely represented as a Laplacian penalization of the objective function. 
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In addition, a Laplacian regularizing term has been used to 
obtain desirable solution characteristics for a variety of other 
optimization problems (e.g., machine learning problems ED, 
EH, sensor network localization problems EH, and analyses 
of flow networks Ell). 

C. An Algorithm for Determining the Laplacian Weights 

Having established a weighted Laplacian form for the 
objective function, we introduce an iterative algorithm for 
determining appropriate weights D for obtaining a solution 
to (foi l with rank(W) = 1. We note that the proposed 
algorithm is similar in spirit to the method in ll35l Section 2.4], 
which iteratively updates weighting parameters to promote 
low-rank solutions of SDPs related to image reconstruction 
problems. 

The proposed algorithm is inspired by the apparent power 
line flow penalty used in CD and the iterative approach 
to determining appropriate buses for enforcing higher-order 
moment constraints in M- The approach in iTT^ penalizes 
the apparent power flows on lines associated with certain 
submatrices of W that are not rank one3 Similar to the 
approach in mi, the proposed algorithm adds terms to the 
objective function that are associated with certain “problematic 
lines.” 

The heuristic for identifying problematic lines is inspired 
by the approach used in 03 to detect “problematic buses” 
for application of higher-order moment constraints. Denote the 
solution to (fTSl i as W* and the closest rank-one matrix to W* 
as (By the Eckart and Young theorem E3, the closest 

rank-one matrix is calculated using the eigendecomposition 
W^t) = Xir/ipJ, where Ai and rji are the largest eigenvalue of 
W* and its associated unit-length eigenvector, respectively.) 
If W* = W(t), then rank(W*) = 1 and we can recover 
the global optimum to El using (nni. Otherwise, previous 
work lfT3 compares the power injections associated with W* 
and to calculate power injection mismatches S'™-’ 

for each bus k G Af: 

mis _ 

tr {Yfe (W* - } + jtr {y^ (w* - } 

(19) 

where | • | denotes the magnitude of the complex argument. In 
the parlance of 03, problematic buses are those with large 
power injection mismatches 

To identify problematic lines rather than buses, we 
modify (HSll to calculate apparent power flow mismatches 
for each line {I, m) G £: 

o flow mis 

S(l,m) = 

tr (W* - } + jtr (w* - } 

+ |tr {Zmi (W* - } +]tv\Zmi (W* - } . 

(20) 

*The submatrices are determined by the maximal cliques of a chordal 
supergraph of the network; see 1101 . 1191 . 1231 for further details. 


Algorithm 1 Iterative Algorithm for Determining Weights 
1: Input; tolerances efiow and Cinj, max relaxation gap 5 
2: Set D = Onixni 

3: Solve the SDP relaxation (|9|l to obtain c* 

4: Calculate and using (|20l) and (O 

5: while termination criteria not satisfied 
6: Update weights; D G- D -|- diag 

7: Solve the generation-cost-constrained relaxation (fT3]) 

8: Calculate and using (|20ll and (UHl 

9: end while 

10: Calculate the voltage phasors using (fTOl) and terminate 


Observe that sums the magnitude of the apparent 

power flow mismatches at both ends of each line. 

The condition rank(W*) = 1 (i.e., “feasibility” in this 
context) is considered satisfied for practical purposes using the 
criterion that the maximum line flow and power injection mis- 
matches (i.e., max(i_^)g£ and maxfcgAf Jfc ) 

are less than specified tolerances efiow and emj, respectively, 
and the voltage magnitude limits (l3dl are satisfied to within a 
specified tolerance ey0 

As described in Algorithm [D the weights on the diagonal of 
D are determined from the line flow mismatches gfjowmis ^ 

(Z,m) 

Specifically, the proposed algorithm first solves the SDP 
relaxation (|9l) to obtain both the lower bound c* on the optimal 
objective value and the initial line flow and power injection 
mismatches V {l,m) G C and vfc G N. 

While the termination criteria (maxp ^ 

efiow, maxfeg^vy < einj, and no voltage limits 

violated by more than ey) are not satisfied, the algorithm 
solves (El (i.e., the SDP relaxation with the constraint ensur¬ 
ing that the generation cost is within 5 of the lower bound). 
The objective function is defined using the weighting matrix 
D = diag where diag (•) denotes the matrix with 

the vector argument on the diagonal and other entries equal 
to zero. Each iteration adds the line flow mismatch vector 
g flow mis fj-Qjjj (}jg solution to (foi l to the previous weights 
(i.e., D G- D + diag 

Upon satisfaction of the termination criteria, the algorithm 
uses ([Tol l to recover a feasible solution to ([D that has an objec¬ 
tive value within <5 of the global optimum. Again, “feasibility” 
in this context is judged using the termination criteria efiow, 
einj, and ey. 

Note that Algorithm [D is not guaranteed to converge. Non¬ 
convergence may be due to the value of S being too small (i.e., 
there does not exist a rank-one solution that satisfies (O) or 
failure to And a rank-one solution that does exist. To address 
the former case. Algorithm [D could be modified to include an 
“outer loop” that increments d by a specified amount (e.g., 
0.5%) if convergence is not achieved in a certain number of 
iterations. We note that, like other convex relaxation methods, 
the proposed approach would benefit from further theoretical 
work regarding the development of a priori guarantees on the 
size of the relaxation gap for various classes of OPE problems. 

^For all test cases, the voltage magnitude limits were satisfied whenever 
the power injection and line flow mismatch tolerances were achieved. 











7 


For some problems with large relaxation gaps (e.g., the 
3-bus system in 12811 . the 5-bus system in 1^ . and the 9-bus 
system in mx no purely penalization-based methods have so 
far successfully addressed the latter case where the proposed 
algorithm fails to find a rank-one solution that satisfies the 
generation cost constraint (fTSl i with sufficiently large 6 (i.e., 
no known penalty parameters yield feasible solutions using 
the methods in mi, mi for these test cases). One possible 
approach for addressing this latter case is the combination 
of penalization techniques with Lasserre’s moment relaxation 
hierarchy HD-IH. The combination of the moment relax¬ 
ations with the penalization methods enables the computation 
of near-globally-optimal solutions for a broader class of OFF 
problems than either method achieves individually. See CSl 
for further details on this approach. 

We note that despite the lack of a convergence guarantee, the 
examples in Section|V]demonstrate that Algorithm[T]is capable 
of finding feasible points that are near the global optimum for 
many OFF problems, including large test cases. 

V. Results 

This section demonstrates the effectiveness of the proposed 
approach using several small example problems as well as 
large test cases representing portions of European power 
systems. The SDF relaxation yields a small but non-zero 
relaxation gap for the test cases selected in this section, and 
Algorithm [T] yields points that are feasible for Q (to within 
the specified termination criteria) and that are near the global 
optimum for these test cases. For other test cases with a large 
SDF relaxation gap, such as those mentioned earlier in m, 
Ea-®, the proposed algorithm does not converge when 
tested with a variety of values for S. 

The results in this section use line flow and power injection 
mismatch tolerances efiow and einj that are both equal to 
1 MVA and ey = 5 x 10“"^ per unit. The implementation of 
Algorithm[I]uses MATLAB 2013a, YALMIF 2015.02.04 ISl, 
and Mosek 7.1.0.28 138], and was solved using a computer 
with a quad-core 2.70 GHz processor and 16 GB of RAM. 

Applying Algorithm [T] to several small- to medium-size 
test cases from ifT^ . ll20l . ll26l . ll2^ . yields the results 
shown in Table |I] Tables |II] and |III] show the results from 
applying Algorithm [T] to large test cases which minimize 
generation cost and active power losses, respectively. These 
test cases, which are from m and ll25l . represent portions 
of European power systems. The SDF relaxation (|9]l has a 
small but non-zero relaxation gap for all test cases consid¬ 
ered in this section. The columns of Tables HI-IIIII show the 
case name and reference, the number of iterations of Algo¬ 
rithm [T] the final maximum apparent power flow mismatch 
maxp MVA, the final maximum power 

injection mismatch } in MVA, the specihed 

value of 6, an upper bound on the relaxation gap from the 
solution to the SDF relaxation and the total solver time 
in seconds. 

Note that the large test cases in Tables HI] and |In| were 
preprocessed to remove low-impedance lines as described 
in El in order to improve the numerical convergence of 


TABLE I 

Results for Small and Medium Size Test Cases 


Case 

Name 

Num. 

Iter. 

Max 

Flow Mis. 
(MVA) 

Max 
Inj. Mis. 

(MVA) 

<5 

(%) 

Max 

Relax. 

Gap (%) 

Solver 

Time 

(sec) 

LMBD3 1201 

1 

1.3E-5 

1.6E-5 

0.5 

0.50 

0.7 

MLD.3 1291 

1 

7.3E-6 

7.2E-5 

0.5 

0.50 

0.5 

MHI4Q IHI 

2 

1.8E-5 

9.9E-6 

0.5 

0.02 

1.2 

MHI4L II4l 

2 

8.1E-5 

7.8E-5 

0.5 

0.33 

1.2 

KDS14Lin |30| 

1 

1.2E-3 

9.2E-4 

1.0 

1.00 

0.7 

KDSMQuad |3()| 

1 

1.4E-4 

8.4E-5 

1.0 

1.00 

0.6 

KDS30Lin |30| 

7 

9.3E-1 

9.2E-1 

2.5 

2.50 

4.6 

KDS30Quad |30| 

6 

8.1E-1 

8.0E-1 

2.0 

2.00 

3.6 

KDS30IEEEOuad 1301 

100 

9.5E-1 

7.2E-1 

2.5 

2.50 

129.8 

MH39L II4l 

1 

1.3E-2 

9.8E-3 

0.5 

0.27 

0.7 

MH570 Il4l 

1 

1.2E-3 

6.9E-4 

0.5 

0.03 

0.7 

MH57L II4l 

1 

3.2E-4 

5.2E-4 

0.5 

0.16 

0.9 

MHII80 1171 

2 

3.3E-3 

2.7E-3 

0.5 

0.50 

2.6 

MHII8L 11^ 

2 

3.1E-3 

3.1E-3 

1.0 

1.00 

3.3 

IEEE 300 |2^ 

1 

1.3E-1 

1.2E-1 

0.5 

0.01 

3.0 


the SDF relaxation. Lines which have impedance magnitudes 
less than a threshold (thrshz in ifTSl l of 1 X 10 ^ per unit 
are eliminated by merging the terminal buses. Table HYI 
describes the number of buses and lines before and after 
this preprocessing. Low-impedance line preprocessing was not 
needed for the test cases in Table U After preprocessing, 
MOSEK’s SDF solver converged with sufficient accuracy to 
yield solutions that satisfied the voltage magnitude limits to 
within ey = 1 X 10“^ per unit and the power injection and 
line flow constraints to within the corresponding mismatches 
shown in Tables aHmi 

These results show that Algorithm [T] finds feasible points 
(within the specified tolerances) that have objective values 
near the global optimum for a variety of test cases. Further, 
Algorithm [T] globally solves all OFF problems for which the 
SDF relaxation (|9]) is exact (e.g., many of the IEEE test 
cases m, several of the Folish test systems m, and the 89- 
bus FEGASE system li25l ). Thus, the algorithm is a practical 
approach for addressing a broad class of OFF problems. 

We note, however, that Algorithm [T] does not yield a 
feasible point for all OFF problems. For instance, the test 
case WB39mod from m has line flow and power injection 
mismatches of 18.22 MVA and 12.99 MVA, respectively, after 
1000 iterations of Algorithm [T] The challenge associated with 
this case seems to result from light loading with limited ability 
to absorb a surplus of reactive power injections, yielding 
at least two local solutions. In addition to challenging the 
method proposed in this paper, no known penalty parameters 
yield feasible solutions to this problem. Generalizations of the 
SDF relaxation using the Lasserre hierarchy have successfully 
calculated the global solution to this case m, Us]. Further, 
while Algorithm [T] converges for five of the seven test cases 
in ll^ which have small relaxation gaps (less than 2.5%), 
the algorithm fails for two other such test cases as well as 
several other test cases in ll^ which have large relaxation 
gaps. We note that the tree topologies used in the test cases 
in lIMl are a significant departure from the mesh networks used 
in the standard test cases from which they were derived; the 
proposed algorithm succeeds for several test cases that share 
the original network topologies. 

Note that for the large test cases in Tables [Hi and |ni| 
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TABLE II 

Results for Large Test Cases that Minimize Generation Cost 


Case 

Name 

Num. 

Iter. 

Max 

Flow Mis. 
(MVA) 

Max 
Inj. Mis. 

(MVA) 

S 

(%) 

Max 
Relax. 
Gap (%) 

Solver 

Time 

(see) 

PL-2383wp (241 

2 

0.54 

0.50 

0.5 

0.50 

78.6 

PL-30 12wd 1241 

2 

0.36 

0.27 

0.5 

0.50 

107.6 

PL-3120SP 1241 

2 

0.56 

0.33 

0.5 

0.50 

84.2 


TABLE III 

Results for Large Test Cases that Minimize Active Power Loss 


Case 

Name 

Num. 

Iter. 

Max 

Flow Mis. 
(MVA) 

Max 
Inj. Mis. 

(MVA) 

S 

(%) 

Max 

Relax. 

Gap (%) 

Solver 

Time 

(sec) 

PL-2383wp [24J 

5 

0.21 

0.16 

0.5 

0.26 

154.0 

PL-30 12wd 1241 

5 

0.08 

0.04 

0.5 

0.18 

232.2 

PL-3120SD 1241 

5 

0.25 

0.19 

0.5 

0.38 

232.6 

PEGASE-1354 1251 

12 

0.27 

0.18 

0.5 

0.15 

199.2 

PEGASE-2869 1251 

38 

0.91 

0.69 

0.5 

0.15 

2378.4 


TABLE IV 

Descriptions of Large Test Cases Before and After 
Low-Impedance Line Preprocessing 


Case 

Name 

Before Preprocessing 

After Preprocessing 

Num. 

Buses 

Num. 

Lines 

Num. 

Buses 

Num. 

Lines 

PL-2383wp 

2,383 

2,869 

2,177 

2,690 

PL-3012wp 

3,012 

3,572 

2,292 

2,851 

PL-3120sp 

3,120 

3,693 

2,314 

2,886 

PEGASE-1354 

1,354 

1,991 

1,179 

1,803 

PEGASE-2869 

2,869 

4,582 

2,120 

4,164 


Algorithm [T] is often computationally faster and has a 
more straightforward computational implementation than the 
moment-relaxation-based approaches in lfT4ll . lITSl . However, 
Algorithm [T] results in feasible points with larger objective 
values and does not solve as broad a class of OPF problems 
as existing moment-relaxation-based approaches in lfT4l . ITSl. 

Numerical experience suggests that S = 0.5% is usually an 
appropriate parameter choice: as discussed in Section IIV-AI 
the SDP relaxation gap is smaller than 0.5% for many test 
cases. For OPF problems with a significantly larger relaxation 
gap, the proposed approach typically fails to yield a feasible 
solution. Thus, values of 6 that differ significantly from 0.5% 
are not likely to be useful in practice. 

We note that the interior point solver in Matpower ob¬ 
tained superior relaxation gaps for the test cases considered 
in this paper. Within approximately five seconds for the large 
test cases in Tables HJ and HIH Matpower obtained relaxation 
gaps that ranged from 0.14% to 0.32% smaller than those 
obtained with Algorithm [T] (Of course, Matpower cannot 
provide any measure of the quality of its solution in terms of 
a lower bound on the globally optimal objective value whereas 
Algorithm[T]provides such guarantees.) The smaller relaxation 
gaps obtained using Matpower suggest that smaller values of 
i5 could be used in Algorithm [T] Indeed, additional numerical 
experiments demonstrated that Algorithm [T] converged with 
S = 0.25% (half the value used in previous numerical exper¬ 
iments) for all test cases for which the Matpower solution 
indicated that a value of i5 = 0.25% was achievable. 

We select termination parameter values of e flow and einj of 
1 MVA, which is a reasonable value for practical power system 



Eig. 1. Maximum Apparent Power Flow Mismatches versus Iteration of 
Algorithm [T] for Active Power Loss Minimizing Test Cases 


Maximum Active and Reactive Fiow Mismatch 



Pig. 2. Maximum Active and Reactive Power Plow Mismatches versus 
Iteration of Algorithm [T] for Active Power Loss Minimizing Test Cases 

applications. This tolerance is typically numerically achievable 
with MOSEK’s SDP solver, which experience suggests is often 
a limiting factor to obtaining smaller mismatches. 

Note that the maximum mismatches do not necessarily 
decrease monotonically with each iteration of Algorithm [T] 
Figs. [I] and |2] show the maximum flow mismatches (on a 
logarithmic scale) for the test cases that minimize active 
power losses (cf. Table HIH i. Likewise, Figs. |3] and |4] show 
the maximum power injection mismatches for the same test 
cases. Although the mismatches do not always decrease mono¬ 
tonically, there is a generally decreasing trend which results 
in satisfaction of the termination criteria for each test case. 
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Fig. 3. Maximum Apparent Power Injection Mismatches versus Iteration of 
Algorithm ^ for Active Power Loss Minimizing Test Cases 



Fig. 4. Maximum Active and Reactive Power Injection Mismatches versus 
Iteration of Algorithm ^ for Active Power Loss Minimizing Test Cases 

At each iteration. Algorithm [T] yields larger reactive power 
mismatches than active power mismatches for these test cases. 

Note that it is not straightforward to compare the compu¬ 
tational costs of the Laplacian objective approach and other 
penalization approaches in the literature flSl . jT^ . A single 
solution of the penalized SDP relaxations in CD requires 
approximately the same computational effort as one iteration 
of Algorithm [T] Thus, if one knows appropriate penalty 
parameters, the method in lfT9l is faster for problems where 
the SDP relaxation is not exact. However, the key advantage of 
the proposed approach is that there is no need to specify any 
parameters other than the value of S used in the generation cost 


constraint. In contrast, the literature largely lacks systematic 
approaches for identifying appropriate parameter values for 
the penalization methods in ifTSl , ifT^ . 

VI. Conclusion 

The SDP relaxation of lO is capable of globally solving 
a variety of OPF problems. To address a broader class of 
OPF problems (i.e., problems for which the SDP relaxation 
is not exact but provides lower bounds that are close to the 
global optima), this paper has described an approach that 
finds feasible points with objective values that are within 
a specified percentage of the global optimum. Specifically, 
the approach in this paper adds a constraint to ensure that 
the generation cost is within a small specified percentage 
of the lower bound obtained from the SDP relaxation. This 
constraint frees the objective function to be chosen to yield 
a feasible (i.e., rank-one) solution rather than a minimum- 
cost solution. Inspired by previous penalization approaches 
and results in the optimization literature, an objective function 
based on a weighted Laplacian matrix is selected. The weights 
for this matrix are iteratively determined using “line flow 
mismatches.” The proposed approach is validated through 
successful application to a variety of both small and large 
test cases, including several OPF problems representing large 
portions of European power systems. There are, however, test 
cases for which the approach takes many iterations to converge 
or does not converge at all. 

Future work includes modifying the algorithm for choos¬ 
ing the weights in order to more consistently require fewer 
iterations. Also, future work includes testing alternative SDP 
solution approaches with “hot start” capabilities to improve 
computational efficiency by leveraging knowledge of the so¬ 
lution to a “nearby” problem from the previous iteration of 
the algorithm. Future work also includes extension of the 
algorithm to a broader class of OPF problems, such as the 
test case WB39mod from III and several examples in ll30l . 

Additional future work includes leveraging recent re¬ 
sults showing that constraints from alternative relaxations 
(e.g., HD-lUl, IE), EQI, ED, EOl) can tighten the SDP 
relaxation. Augmenting the proposed approach with such 
constraints may increase its applicability to an broader class 
of problems. 
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